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Abstract — In cognitive radio systems, secondary users can 
utilize multiple dispersed bands that are not used by primary 
users. In this paper, time delay estimation of signals that occupy 
multiple dispersed bands is studied. First, theoretical limits on 
time delay estimation are reviewed. Then, two-step time delay 
estimators that provide trade-offs between computational com- 
plexity and performance are investigated. In addition, asymptotic 
optimality properties of the two-step time delay estimators are 
discussed. Finally, simulation results are presented to explain the 
theoretical results. 

Index Terms — Cognitive radio, time delay estimation, 
maximum-likelihood, diversity, Cramer-Rao lower bound. 

I. Introduction 

Cognitive radio presents a promising approach to implement 
intelligent wireless communications systems 0~)-[|4]. Cognitive 
radios can be regarded as more capable versions of software 
defined radios in the sense that they have sensing, awareness, 
learning, adaptation, goal driven autonomous operation and 
reconfigurability features |5), |6), which facilitate efficient 
use of radio resources such as power and bandwidth QJ. 
As the electromagnetic spectrum is a precious resource, it is 
important not to waste it. The recent spectrum measurement 
campaigns in the United States Q and Europe [8 1 indicate 
that the spectrum is under-utilized; hence, opportunistic use 
of unoccupied frequency bands is highly desirable. 

Cognitive radio presents a solution to inefficient spectrum 
utilization by opportunistically using the available spectrum of 
a legacy system without interfering with the primary users of 
that spectrum J2), J3). In order to facilitate such opportunistic 
spectrum utilization, cognitive radio devices should be aware 
of their locations, and monitor the environment continuously. 
Location awareness requires a cognitive radio device to per- 
form accurate estimation of its position. Cognitive radio de- 
vices can obtain position information based on the estimation 
of position related parameters of signals traveling between 
them |9), iflOl . Among various position related parameters, the 
time delay parameter provides accurate position information 
with reasonable complexity IflOl . The main focus of this paper 
is time delay estimation in cognitive radio systems. 

The main difference between time delay estimation in 
cognitive radio systems and that in conventional systems is that 
a cognitive radio system can transmit and receive over multiple 
dispersed bands. In other words, as a cognitive radio device 
can use the spectral holes in a legacy system, it can have a 
spectrum that consists of multiple bands that are dispersed over 
a wide range of frequencies (cf. Fig. [TJ. In QT), the Cramer- 
Rao lower bounds (CRLBs) for time delay estimation are 
obtained for dispersed spectrum cognitive radio systems, and 
the effects of carrier frequency offset (CFO) and modulation 
schemes of training signals on the accuracy of time delay 
estimation are quantified. The CRLB expressions imply that 
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frequency diversity can be utilized in time delay estimation. 
Similarly, the effects of spatial diversity on time delay estima- 
tion are investigated in ITZll for single-input multiple-output 
systems. Also, the effects of multiple antennas on time delay 
estimation and synchronization problems are studied in lfl3l . 

This paper studies time delay estimation for dispersed 
spectrum cognitive radio systems. First, the theoretical limits 
on time delay estimation are reviewed, and the concept of fre- 
quency diversity for time delay estimation is discussed. Then, 
optimal and suboptimal time delay estimation techniques are 
studied. Since optimal maximum likelihood (ML) time delay 
estimation can have very high computational complexity for 
signals with multiple dispersed bands, two-step time delay 
estimation techniques are investigated. The proposed two-step 
time delay estimators first extract unknown parameters related 
to signals in different frequency bands, and then obtain the 
final time delay estimate in the second step. In other words, 
multiple observations (signals at different frequency bands) 
are processed efficiently to provide a trade-off between com- 
putational complexity and estimation performance. In addition, 
the optimality properties of the two-step estimators are inves- 
tigated for high signal-to-noise ratios (SNRs), and simulation 
results are presented to verify the theoretical analysis. 

II. Signal Model 

Consider a scenario in which K dispersed frequency bands 
are available to the cognitive radio system, as shown in Fig. 
Q] The transmitter generates a signal that occupies all the K 
bands simultaneously, and sends it to the receiver. Then, the 
receiver is to estimate the time delay of the incoming signal. 
Since the available bands can be quite dispersed, the use of 
orthogonal frequency division multiplexing (OFDM) approach 
lfl4l can require processing of very large bandwidths. There- 
fore, processing of the received signal in multiple branches is 
considered in this study, as in Fig. [2] ifTTl . 

For the receiver in Fig. |2j the baseband representation of 
the received signal in the ith branch is given by 



n(t) = a i e J ' Uit s i (i-T) +rn{t) , 



(1) 



for i 



1. 



, K, where r is the time delay of the signal, 
a, = a, c?'^ i and u>i represent, respectively, the channel 
coefficient and the CFO for the signal in the ith branch, Si(t) 
is the baseband representation of the transmitted signal in the 
ith band, and n.j(i) is complex white Gaussian noise with 
independent components, each having spectral density of . It 
is assumed that the signal in each branch can be modeled 
as a narrowband signal; hence, a single complex channel 
coefficient is used to represent the fading of each signal. Also, 
it should be noted that the effects of CFO are considered in 
the signal model in (HJ since multiple down-conversion units 
are employed in the receiver, as shown in Fig. |2] 
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Fig. 1 

Illustration of dispersed spectrum utilization in the cognitive 
radio system, where the white spaces represent the available 

BANDS. 




Jc K 

Fig. 2 

Block diagram of the front-end of the cognitive radio 
receiver. 



III. Theoretical Limits on Time Delay Estimation 

Estimation of the time delay parameter r based on K 
received signals in (Q3 involves (3K+ 1) unknown parameters 
since the channel coefficients and CFOs are also unknown. 
In other words, the vector of unknown parameters, 9, can be 
expressed as 9 = [t a% ■ ■ ■ o,k <fri ■ ■ ■ <Pk wi • • ■ uk] ■ 

For an observation interval of [0, T], the log-likelihood 
function for 9 is expressed a£] |[T5l 

A(9) =c-J2t^ Oi^Siit -t)\ 2 dt, (2) 

where c is a constant that is independent of 9. Then, the Fisher 
information matrix (FIM) lfT31 can be obtained from (f2]i as 
in ifTTl . and the inverse of the FIM can be used to obtain 
the CRLB on mean-squared errors (MSEs) of unbiased time 
delay estimators. In its most generic form, the CRLB can be 
expressed as IfTTl 

E{(f - r) 2 } > (E ^ - (4 R )7^) - tj , O) 

where Ei = J Q T \si(t — r)\ 2 dt is the signal energy, Ei = 
Jo l s '(^ — r )| 2al i, with s'(t) representing the first derivative 
of s(t), and Ef = ' K{s' % (t - n)s*{t - n)}dt, with TZ 
denoting the operator that selects the real-part of its argument. 
In addition, £ represents a term that depends on the spectral 
properties of signals Sj(t) for i = 1, . . . , K IfTTl . 

The CRLB expression in (O reveals that the accuracy of 
time delay estimation depends on the SNR at each branch (via 
the aj /of terms), as well as on the properties of signals Sj (t) 
in ([Tj. In addition, the summation term in (0 indicates that 
the accuracy can be improved as more bands are employed, 
which implies that frequency diversity can be utilized in time 
delay estimation. For instance, when one of the bands is in a 
deep fade (that is, small af ), some other bands can still be in 

'The unknown parameters are assumed to be constant for t G [0,T]. 
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Fig. 3 

The block diagram of the two-step time delay estimation 
approach. The signals ri(t), . . . , r K (t) are as shown in Fig. [2] 



good condition to facilitate accurate time delay estimation. 

In order to investigate the effects of signal design on the 
time delay estimation accuracy, suppose that the baseband 
representation of the signals in different branches are of the 
form Sj(t) = J2idi,iPi(t — ITi), where du denotes the 
complex training data and is a pulse with duration Tj. 
Then, the £ term in (0) becomes zero, which results in a CRLB 
expression that would be obtained in the absence of CFOs 
fTJJ. In other words, the effects of CFOs can be mitigated via 
appropriate signal design. 

In the special case of \dn\ = V7 and pAt) satisfying 
Pi (0) = pi(Ti) for % = 1, . . .', K, © reduces to fill 

E{(f-r) 2 }>(£^#) . (4 ) 

Hence, for linearly modulated signals with constant envelopes, 
improved time delay estimation accuracy can be achieved. 

IV. Two-Step Time Delay Estimation 
The ML estimate of 9 can be obtained from (ffji as 

0Mh = arg max / TZ I a*e~ 3Uit ri(t)sl(t - r) \ dt - L 

9 ~i a i Jo 1 ' 2cr, 

(5) 

which requires an optimization over a (3K + l)-dimensional 
space, hence is quite impractical in general. Therefore, a two- 
step time delay estimation approach is considered in this study, 
as shown in Fig. [3] In the first step, each branch of the receiver 
performs estimation of the time delay, the channel coefficient 
and the CFO related to the signal in that branch. Then, in the 
second step, the estimates from all the branches are used to 
obtain the final time delay estimate. 

A. First Step: Parameter Estimation at Different Branches 

In the first step, the unknown parameters of each received 
signal are estimated at the corresponding receiver branch 
according to the ML criterion (cf. Fig.0. Based on the signal 
model in ([JJ, the log-likelihood function at branch i becomes 

1 f T 

A<(0<) = <k - —j / \n(t) - at &> Uit Si(t - r)| 2 dt , (6) 
l0 ~i Jo 

for i = 1, . . . , K, where 9i = [r cij 4>i uji] denotes the vector 
of unknown parameters related to the signal at the ith branch, 
Ti(t), and a is a constant that is independent of 9i. 
From ©, the ML estimator at branch i is expressed as 

9i = argiriin / k(t) - a t e J ^ t s l {t - r)| 2 dt , (7) 

°i JO 



where 6i = [fi hi 4>i u>i\ is the vector of estimates at the ith 
branch. After some manipulation, (|7]i yields 

K{r i {t)e-X u < 1 +M8' i (t-Ti)}dt 

(8) 



arg max 
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{n(t)e-^* t+ ^8t(t-n)}dt . 



(9) 



In other words, at each branch, optimization over a three- 
dimensional space is performed to obtain the unknown pa- 
rameters, which is significantly less complex than the ML 
estimation in (0 that requires optimization over (3K + 1) 
variables. 

B. Second Step: Combining Estimates from Different Branches 

After obtaining time delay estimates fi, . . . ,tk in ®, the 
second step combines those estimates according to one of the 
criteria below and makes the final time delay estimate. 

1 ) Optimal Combining: According to the "optimal" com- 
bining] criterion, the time delay estimate is obtained as 



r = 



(10) 



where f,; is the time delay estimate of the ith branch, which is 
obtained from ([8]), and Ki = afEi/af. In other words, the 
optimal combining approach estimates the time delay as a 
weighted average of the time delays at different branches, 
where the weights are chosen as proportional to the mul- 
tiplication of the SNR estimate, Eid,?/af, and Ei/Ei. As 
Ei is defined as the energy of the first derivative of Si(t), 
Ei/Ei can be expressed, using Parseval's relation, as Ei/Ei = 
47r 2 /3 2 , where Pi is the effective bandwidth of s,(4), which is 
defined as (3f = ± / 2 |^(/)| 2 d/ with S t (f) denoting 
the Fourier transform of Si(t) lfl5l . Therefore, it is observed 
that the optimal combining technique assigns a weight to 
the time delay estimate of a given branch in proportion to 
the product of the SNR estimate and the effective bandwidth 
related to that branch. The intuition behind this combining 
approach is that signals with larger effective bandwidths and/or 
larger SNRs facilitate more accurate time delay estimation 
JT5J ; hence, their weights are larger in th e com bining process. 
This intuition will be verified in Section IIV-CI theoretically. 

2) Selection Combining (SC): Another approach to obtain 
the final time delay estimate is to determine the "best" branch 
and to use its estimate as the final time delay estimate. 
According to SC, the best branch is defined as the one with 
the maximum value of Ki = af Ei/af for i — 1, . . . , K. That 
is, the branch with the maximum multiplication of the SNR 
estimate and the effective bandwidth is selected as the best 
branch and its estimate is used as the final one. In other words, 



ie{i,...,K} I J 



1) 



where f rn represents the time delay estimate at the mth branch. 

3) Equal Combining: The equal combining approach as- 
signs equal weights to the estimates from different branches 
and obtains the time delay estimate as r = j? Yld=i ^ ■ 



Considering the combining techniques above, it is observed 
that they are counterparts of the diversity combining tech- 
niques employed in communications systems |[T6l . However, 
the main distinction is that the aim is to maximize the SNR or 
to reduce the probability of symbol error in communications 
systems |[T6l . whereas, in the current problem, it is to reduce 
the MSE of the time delay estimation. In other words, this 
study focuses on diversity combining for time delay estima- 
tion, where the diversity is due to the dispersed spectrum 
utilization of the cognitive radio system. 

C. Optimality Properties of Two-Step Time Delay Estimation 

In this section, the asymptotic optimality properties of the 
two-step time delay estimators are investigated in the absence 
of CFO. In order to analyze the performance of the estimators 
at high SNRs, the result in flTZI for time-delay estimation at 
multiple receive antennas is considered first. 

Lemma 1 [12): Assume that s-(< — r)s*(t - r)dt = 
for i = 1, . . . , K. Then, for the signal model in (fJJ, the delay 
estimate in ((S]) and the channel amplitude estimate in ((9]) can 
be modeled, at high SNR, as 



+ Vi and &i = a,i + rji 



(12) 



for i = 1, . . . , K, where Vi and rji are independent zero mean 
Gaussian random variables with variances af/(Eiaf) and 
a\/Ei, respectively. In addition, z/j and Vj (rji and rjj) are 
independent for i =/= j. 

From Lemma 1, it is obtained that E{fi} = r for i = 
1,...,K. In other words, the time delay estimates of all 
branches are asymptotically unbiased. Since the combining 
techniques in the previous section considers only one, or a 
linear combinations of the time delay estimates at different 
branches, the two-step time delay estimation techniques have 
an asymptotic unbiasedness property. 

Regarding the variance of the estimators, it is first shown 
that the optimal combining technique has a variance that is 
approximately equal to the CRLB at high SNRsQ To that 
aim, the conditional variance of f in dTOb given fli, . . . , clk 
is expressed as follows: 



Yn=i K i Var{f<|ai, . . . ,a K ] 



Var{f |oi, . . .,a K } 



(13) 



where the independence of the time delay estimates is used to 
obtain the result (cf. Lemma 1). Since Var{fj|di, . . . , clk} = 
Var{f.;|di} = o-f/(Eiaf) from Lemma 1 and Ki = ajEi/a? , 
(fT3l can be manipulated to obtain 



(14) 



Lemma 1 states that at is distributed as a Gaussian random 
variable with mean a, and variance cr^/Ei at high SNRs. 
Hence, for sufficiently large values of , ( TT4l can 



be approximated by ( YliLi ~i£ L ) > which is equal to the 
CRLB expression in ©. Therefore, the optimal combining 
technique in (TTOb yields an approximately optimal estimator 
at high SNRs. 



2 The optimality property is investigated in Section llV-CI 



3 This is the main reason why this combining technique is called optimal. 
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Fig. 4 

RMSE VS. SNR FOR THE TWO-STEP ALGORITHMS, AND THE 
THEORETICAL LIMIT (CRLB). THE SIGNAL OCCUPIES THREE DISPERSED 
BANDS WITH B\ = 200 KHZ, B 2 = 100 KHZ AND B 3 = 400 KHZ. 
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Fig. 5 

RMSE VS. THE NUMBER OF BANDS FOR THE TWO-STEP ALGORITHMS, 
AND THE THEORETICAL LIMIT (CRLB). EACH BAND IS 100 KHZ WIDE, 
AND of = 0.1 Vi. 



Regarding the selection combining approach in ( fTTT i. the 
conditional variance can be approximated at high SNR as 



Var{f|a!, . . . , a^-} w min • 



K 



Bkal 



E K a 2 K 



(15) 



In general, the SC approach performs worse than the optimal 
combining technique. However, when the estimate of a branch 
is significantly more accurate than the others, its performance 
can get very close to that of the optimal combining technique. 

Final ly, for the equal combining technique in Section 
IIV-B.3I the variance can be calculated as Var{f} = 

■^r X^Li ■ The equal combining approach is expected 

to have the worst performance since it does not make use of 
any information about the SNR or the signal bandwidths in 
the estimation of the time delay, as investigated next. 

V. Simulation Results and Conclusions 

In this section, simulations are performed to evaluate the 
CRLBs and the performance of the time delay estimators. 
Signal Si(t) in (Q} at branch i is modeled by a unit-energy 
Gaussian doublet as in |11| with bandwidth B{. In all the 



simulations, the spectral densities of the noise at different 
branches are assumed to be equal; that is, crj = a for i = 
1, . . . , K. Also, the SNR of the system is defined with respect 
to the total energy of the signals at different branches, i.e., 

SNR=101og 10 (£|^). 

In assessing the root-mean-squared errors (RMSEs) of the 
different estimators, a Rayleigh fading channel is considered. 
Namely, the channel coefficient a, = as e- 7 ^' in (UJ is modeled 
as di being a Rayleigh distributed random variable and fa 
being uniformly distributed over [0, 2n). In addition, the same 
average power is assumed for all the bands; that is, E{|a^ | 2 } = 
1 is used. The time delay, r, in ([TJ is uniformly distributed 
over the observation interval, and it is assumed that there is 
no CFO in the system. 

First, the performance of the two-step estimators is eval- 
uated with respect to the SNR for a system with K = 3, 
B x = 200 kHz, B 2 = 100 kHz and B 3 = 400 kHz. The results 
in Fig. @] indicate that the optimal combining technique has the 
best performance as expected from the theoretical analysis, and 
SC, which estimates the delay according to (TTTb . has perfor- 
mance close to that of the optimal combining technique. On 
the other hand, the equal combining technique has significantly 
worse performance than the others, as it combines all the delay 
estimates equally. Since the delay estimates of some branches 
can have very large errors due to fading, the RMSEs of equal 
combining become quite significant. Finally, it is observed 
that the performance of the optimal combining technique gets 
quite close to the CRLB at high SNRs, in agreement with the 
asymptotic arguments in Section HV-CI 

Next, the RMSEs of the two-step estimators are plotted 
against the number of bands in Fig. [5] where each band is 
assumed to have 100 kHz bandwidth. In addition, the spectral 
densities are set to of = a 2 = 0.1 Vi. From Fig. q\ it is 
observed that the optimal combining has better performance 
than the selection combining and equal combining techniques. 
Also, as the number of bands increases, the amount of 
reduction in the RMSE per additional band decreases (i.e., 
diminishing return). In fact, the selection combining technique 
seems to converge to a constant value for large numbers of 
bands. This is intuitive as the selection combining technique 
always uses the estimate from one of the branches; hence, in 
the presence of a sufficiently large number of bands, additional 
bands do not result in a significant increase in the diversity. On 
the other hand, the optimal combining technique has a slope 
that is quite similar to that of the CRLB; that is, it makes use 
of the frequency diversity efficiently. 
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